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The phonon dispersions of Ndi,86Ceo.i4Cu04+i along the 0, 0] direction have been determined 
by inelastic x-ray scattering. Compared to the undoped parent compound, the two highest longitudi- 
nal phonon branches, associated with the Cu-0 bond-stretching and out-of-plane oxygen vibration, 
are shifted to lower energies. Moreover, an anomalous softening of the bond-stretching band is ob- 
served around q = (0.2, 0, 0). These signatures provide evidence for strong electron-phonon coupling 
in this electron-doped high-temperature superconductor. 



While the coupling between electrons and phonons 
is known to be the driving mechanism for Cooper-pair 
formation in conventional superconductors, its role in 
the high-critical-temperature superconductors (HTcS) 
is the subject of intense research efforts. Recently, 

] evidence for electron-phonon coupling has been in- 
voked in the interpretation of inelastic neutron scat- 
tering (INS) and angle-resolved photoemission spec- 
troscopy (ARPES) experiments. The INS studies, car- 

] ried out on Lai.85Sro.i5Cu04+5 0, |[ |L i|, oxygen- 
doped La2Cu04+5 [|j and YBa^CusOg+a l\ reveal an 

I anomalous softening with doping of the highest longitu- 
dinal optical (LO) phonon branch, in particular along the 

', 9 = Oj 0] direction. This branch is assigned to the Cu- 
O bond-stretching mode ||, 0|. The observed softening 

[ has been interpreted as a signature of a strong electron- 
phonon coupling Pi, |3| , which has been discussed since the 

I discovery of HTcS Furthermore, in an energy range 

similar to the LO bond-stretching phonon band, ARPES 

" studies on three different families of hole-doped HTcS 
reveal a distinct "kink" anomaly in the quasiparticle dis- 

I persion [0. The scenario emerging from the above INS 
and ARPES works suggests a strong coupling between 
the charge carriers and the Cu-0 bond-stretching phonon 
modes to be ubiquitous in HTcS materials, at least for 
hole-doped compounds. However, its role in the pair- 
ing mechanism remains completely unclear [Q. At the 
moment, it may not even be excluded that the electron- 
phonon interaction is pair-breaking for the d-wave su- 
perconducting order parameter. Therefore, it appears 
very important to analyze the strength of the phonon 
anomalies in as many cuprate families as possible and to 
compare them with their superconducting properties. In 
this context, the electron-doped cuprates are of central 
importance due to the distinct character of the doped 
charges in this material: Cu 3dx2-y2 (O 2p) for n(p)- 
type cuprates |jll|, leading to a very different electronic 
structure |12|. Since the phonon anomalies are related 



to a coupling between charge fluctuations and phonons, 
charges with different character may induce quite differ- 
ent electron-phonon interactions. 

In this Letter, we present an inelastic x-ray scatter- 
ing (IXS) study of the phonon dispersion in the n-type 
cuprates Ndi.86Ceo.i4Cu04+5 (NCCO). Inelastic x-ray 
scattering can overcome the main limitation of inelas- 
tic neutron scattering, i.e. the need for sufficiently large 
single crystals of high chemical and structural quality 
Lateral x-ray beam sizes of few tens of /im are rou- 
tinely obtained. Moreover, at photon energies around 
10-20 keV and Z > 3, the total cross section is domi- 
nated by photoelectric absorption, and therefore the typ- 
ical x-ray penetration depths for high-Z materials is of 
the order of 10 - 100 /im. Consequently, very small 
samples (down to less than 10^"* mm"^) can be studied 
with signal rates comparable to typical INS experiments 
on cm^^-sized samples. Despite these advantages, little 
work has been done using IXS on the HTcS compounds 
[|l3|. We choose Nd2-xCexCu04+5 for our IXS study, 
since its crystallographic structure is one of the simplest 
among the HTcS, and because extensive INS studies ex- 
ist for its undoped parent compound Nd2Cu04+5 (NCO) 
||l|, ^. Our interest is focused on the [<^, 0,0] direction, 
where the LO branch displays its strongest anomaly for 
hole-doped HTcS The present results reveal that, 

near the zone center, the two highest longitudinal optical 
branches, assigned to the Cu-0 bond-stretching and 0(2) 
vibration modes, are shifted to lower frequencies with re- 
spect to the undoped parent compound. The interpre- 
tation of our data is supported by lattice dynamics cal- 
culations, taking into account a Thomas- Fermi screening 
mechanism. Furthermore, we observe an anomalous soft- 
ening of the highest branch around q = (0.2,0,0). Our 
results demonstrate that this anomalous behavior of the 
high-energy LO phonon branch is a universal property of 
both hole- and electron-doped HTcS compounds, there- 
fore providing further evidence that electron-phonon in- 
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teractions may play an important role in high-Tc super- 
conductivity. Furthermore, the present results are an 
important demonstration of IXS as a powerful tool for 
the study of the lattice dynamics in small, high-quality 
crystals of complex transition metal oxides. 

The experiment was carried out at the very-high- 
energy-resolution IXS beam-line ID16 at the European 
Synchrotron Radiation Facility (ESRF). X-rays from 
an undulator source were monochromated using a Si 
(111) double-crystal monochromator, followed by a high- 
energy-resolution backscattering monochromator [Q, 
operating at 15816 eV (Si (888) reflection order). A 
toroidal gold-coated mirror refocused the x-ray beam 
onto the sample, where a beam size of 250 x 250 /Ltm^ 
full-width-half-maximum (FWHM) was obtained. The 
scattered photons were energy- analyzed by a spherical 
silicon crystal analyzer 3 m in radius, operating at the 
same Bragg reflection as the monochromator [|l5[. The 
total energy resolution was 1.6 THz (6.6 meV) FWHM. 
The momentum transfer Q was selected by rotating the 
3 m spectrometer arm in the scattering plane perpendic- 
ular to the linear x-ray polarization vector of the incident 
beam. The momentum resolution was set to 0.087 
in both the horizontal and the vertical direction by an 
aperture of 20 x 20 mm^ in front of the analyzer. Fur- 
ther experimental details are given elsewhere ([14, |l^ 
and references therein). The sample is a single crystal 
grown by the traveling-solvent floating-zone method in 4 
atm of O2 at Stanford University. It has been reduced 
under pure Ar atmosphere at 920° C for 20 h, followed 
by a further 20 h of exposure at 500°C to a pure O2 
atmosphere. Such a procedure is necessary to produce a 
superconducting phase in Ndi.86Ceo.i4Cu04+5, although 
its exact effect is not understood. Following this treat- 
ment the sample had a narrow superconducting transi- 
tion with an onset temperature of Tc = 24 K. The sample 
is of very good crystalline quality, with a rocking curve 
width of 0.02° (FWHM) around the [/i,0,0] direction. 
It was mounted on the cold finger of a closed-loop he- 
lium cryostat, and cooled to 15 K. The experiment was 
performed in reflection geometry, and the probed scat- 
tering volume corresponded to about 1.5 x 10^"^ mm'^. 
IXS scans were performed in the -2< v <24 THz range, 
in the r = (6, 0, 0) and r = (7, 1, 0) Brillouin zones @. 
The data were collected along the following three lines: 

I) Q = (6 + ^,0,0), in longitudinal configuration (i.e. 
with q=(e, 0,0) HQ), 

II) Q = (7 ± ^, 1, 0) in almost longitudinal configuration 
(z.e. with q= (^,0,0) and (Q • q)/(5 « g), 

III) Q = (7, 1 — ^, 0) in almost transverse configuration 
(q=(0,C,0) and (Q-q)/Q«0). 

The low temperature and high momentum transfer were 
chosen so as to optimize the count rate on the high- 
frequency optical mode while limiting the loss of contrast 
due to the contribution from the tails of the intense low 
frequency acoustic modes. 




FIG. 1: Experimental IXS phonon spectra of 

Ndi.g6Ceo.i4Cu04+i at T = 15 K and corresponding har- 
monic oscillator model best fits (solid and dotted lines). Both 
scans were performed in the r = (7, 1, 0) Brillouin zones with 
a propagation vector of ^ = 0.2 in an almost longitudinal ge- 
ometry along the a* direction (diamonds) and in an almost 
transverse geometry along h* (open circles). 



Fig. ^ shows a typical energy scan in almost longitudi- 
nal geometry at Q = (6.8, 1, 0), corresponding to ^ = 0.2. 
The data are shown together with the results of a fit, 
where the excitations where modeled by harmonic oscil- 
lators, convoluted with the instrumental resolution func- 
tion. Three features can be clearly distinguished: (i) an 
elastic peak at THz, due to chemical disorder and, pos- 
sibly, strain due to the different thermal expansion coeffi- 
cients between the sample and the sample support during 
cooling, (ii) the longitudinal acoustic phonon, centered 
near 2.7 THz and (iii) a weaker feature around 13 THz. 
In the intermediate energy region (from about 4 to 10 
THz), no distinct phonon peaks are resolved due to the 
dominating contribution from the tails of the elastic and 
acoustic phonon signals. The inset of Fig. 1 emphasizes 
the high-energy region of the spectrum. One can clearly 
distinguish two phonons, centered around 12 and 13.5 
THz, respectively. The weak shoulder at around 16 THz 
can be attributed to an admixture of the transverse op- 
tical (TO) mode, as indicated by the comparison (after 
normalization to the same intensity at 20 THz) with the 
equivalent scan in transverse geometry Q = (7,0.8,0). 
Consequently, the two stronger peaks are unambiguously 
assigned to longitudinal optical (LO) modes. From the 
absence of other higher frequency modes up to 23 THz, 
we conclude that the two modes at 12 and 13.5 THz can 
be identified with the 0(2) vibration and Cu-0 bond- 
stretching modes, respectively (see Refs. ||, Q). 

In order to determine the dispersion of the two highest 
LO branches, IXS spectra were recorded for < ^ < 1. In 
Fig. H we show three spectra taken along the (7 — ^, 1, 0) 
direction. Close to the zone center, at q = (0.1,0,0), 
the highest frequency mode is observed slightly above 15 
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FIG. 2: IXS spectra in the r — (7, 1, 0) Brillouin zone with 
propagation vector q along a* , as indicated in the figure. The 
experimental data (circles) are shown together with their cor- 
responding harmonic oscillator model best fits (solid lines), 
as discussed in the text. 



THz. At q = (0.2,0,0) the highest mode is found at 
the much lower energy of 13.5 THz. Finally, for q = 
(0.4,0,0), we again find a high-frequency mode around 
15.5 THz. 

The peak positions extracted from these and many 
other scans are summarized in Fig. |^. The highest 
branch exhibits a sharp dip around q — (0.2,0,0) and 
recovers for larger g-values. This behavior is most likely 
due to an anti-crossing with the second highest branch 
which is mainly associated with vibrations of the 0(2) 
position in the ^-direction. Therefore, within a standard 
anti-crossing framework, one should interpret the high- 
est longitudinal intensities for q = (0.3,0,0) and above 
as being mainly due to 0(2) vibration. Within that sce- 
nario, this second-highest branch increases its frequency 
in the middle of the zone as in the undoped compound, 
and, except for the fact that the Lyddane- Sachs- Teller 
(LO-TO) gap closes, seems to be insensitive to doping. 
The LO bond-stretching mode just above q — (0.2,0,0) 
is then found at quite low energies, ~ 12 THz, but can 
not be unambiguously followed to larger q- values. Never- 
theless, our data document that the LO bond-stretching 
branch in NCCO is strongly renormalized compared to 
undoped NCO, in particular it bends down anomalously 
from the zone center to q = (0.2, 0, 0). 

In order to further validate the correctness of our as- 
signments, we performed a lattice dynamical calculation 
[ p6[ based on a shell model. We used a common poten- 
tial model for cuprates, in which the interatomic poten- 
tials have been derived from a comparison of INS results 
for different HTcS compounds by Chaplot et al. [ p7[ , 
using a screened Coulomb potential, in order to simu- 
late the effect of the free carriers introduced by doping. 
Following Ref. for metallic Lai.85Sro.i5Cu04+5 and 
YBa2Cu307, we replaced the long-range Coulomb po- 
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FIG. 3: Right hand side: (o, •) experimental longitudinal 
phonon frequencies determined from IXS spectra in NCCO 
at T = 15 K along the [5, 0, 0] direction. Solid circles (•) 
emphasize the highest energy frequencies measured. Solid 
(dot-dashed) lines indicate lattice dynamics calculation with 
a screened (unscreened) Coulomb interaction. 
Left hand side: magnification of the high energy portion of 
the right hand side graph showing the dispersion of the top 
two phonon branches. The highest frequency branch (A), 
as measured by INS (from Ref. 0) in the insulating parent 
compound Nd2Cu04+5 is shown for comparison. The exper- 
imental data are in agreement with calculations, except for 
the anomalous softening at g = (0.2, 0, 0) of the two higher 
energy branches (see text). 



tential Vcisi) by K:((j') /£((?), and for the dielectric func- 
tion we take the semi-classical Thomas-Fermi limit of 
e{q) = 1 + K^/q^, where nj. indicates the screening vec- 
tor. The results of the calculation (without screening: 
dot-dashed lines; with screening: solid lines) are shown as 
well in Fig. ^. The lattice dynamics calculations without 
screening have been included, since they reproduce very 
well the experimental dispersion of the undoped parent 
compound 0. The shift at the zone center of the high- 
energy phonon branches of NCCO with respect to NCO 
is due to the closing of a large LO-TO splitting. Indeed, 
the corresponding Ai and A3 branches in Nd2Cu04+5 
are separated by almost 3 THz at the zone center, as ob- 
served by INS . We point out that in our case a strong 
softening due to Thomas-Fermi screening does not imply 
a higher Thomas-Fermi parameter Ks'. actually, we find 
a screening vector of about 0.39 A , which is com- 
parable to that for Lai.85Sro.i5Cu04+5 ijl^. Though our 
modified calculation reproduces the closing of the LO- 



4 



to 



c 

B 



1 1 1 1 1 

a) 


1 ' 1 ' 1 ' 

b) 


- IL 




/ * \ 






1 



0.04 



0.02 



0.25 0.5 



0.75 0.25 

[4.0,0] 



0.5 0.75 1 



FIG. 4: Comparison between experimental (diamonds) and 
calculated (line) phonon intensities for the two highest opti- 
cal phonon branches along (7 — 1,0). a) bond-stretching 
LO mode starting around 15 THz and b) 0(2) vibration LO 
branch starting around 10.5 THz. 



TO splitting, we still observe an anomalous additional 
softening of the highest bond-stretching LO branch near 
q = (0.2,0,0), which is not reproduced by our calcu- 
lations (see Fig. ^). This branch softens in frequency 
from q = (0.1,0,0) to q = (0.2,0,0) by about Ai/ w 
1.5 THz, which is a shift comparable to the anomalous 
shift observed in Lai.85Sro.i5Cu04+5 at slightly larger q 
[0, ^, Therefore, we believe that this anomalous soft- 
ening is of the same nature as the one observed in p-type 
Lai.85Sro.i5Cu04+5. 

A comparison between the experimental and calculated 
intensities for the two highest phonon branches is shown 
in Fig. |[ The good agreement of the observed inte- 
grated intensities with the calculated ones for the upper 
branches validates the correctness of our phonon branch 
assignment, at least for ^ < 0.2. For ^ > 0.2 we would 
have expected an intensity exchange between the two 
highest branches, which seems to be not observed. 

The main difference between Lai.85Sro.i5Cu04+5- and 
Ndi.g6Ceo.i4Cu04+i is, besides the screening effect, that 
in NCCO the Cu-0 bond-stretching branch is closer in 
energy to the out-of-plane oxygen vibration one, hav- 
ing almost the same energy at ^ = 0.2. These two 
branches belong to the same symmetry and therefore 
cannot cross, so that for ^ > 0.2 softening implies in- 
teraction with the out-of-plane oxygen vibration mode 
with the same symmetry. In the region between q — 
(0.25,0,0) and (0.3,0,0), the two modes are poorly de- 
fined in energy, which is consistent with what is observed 
in Lai.85Sro.i5Cu04+5 by Pintschovius and Braden Q| 
and McQueeney et al. for ^ ^ 0.25 — 0.3. The corre- 
sponding real space periodicity of 3 to 4 unit cells may 
therefore be linked to the proximity to some charge in- 
stability. We remark that a reduced vector ^ ^ 0.25 — 0.3 
approximately corresponds to the nesting vector along 
[^00] direction, as can be inferred from the ARPES data 



of Ref. 

In conclusion, the present results reveal that the 
anomalous softening previously observed in hole-doped 
compounds ||l|, |2[ ||, ^, ||, || , is also present in the electron- 
doped cuprates. This is evidenced by the comparison of 
the present results on doped NCCO with the previously 
reported ones on pure NCO [Q. This implies that: (i) the 
anomaly also exists in n-type cuprates, giving strength to 
the hypothesis [0, ||, |[ ^, ||, |^ of an clcctron-phonon cou- 
pling origin of this feature; (ii) this is a generic feature 
of the high-temperature superconductors, as expected, 
if phonons are relevant to high temperature supercon- 
ductivity. Moreover, this Letter demonstrates that high- 
energy resolution inelastic x-ray scattering has developed 
into an invaluable tool for the study of the lattice dynam- 
ics of complex transition metal oxides, allowing measure- 
ments on small high-quality single crystals which are in- 
accessible to the traditional method of inelastic neutron 
scattering. 
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